Energy-efficient motion camouflage in three dimensions 
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Abstract 

Recent observations suggest that one insect may camouflage its own motion whilst tracking another. 
Here we present a strategy by which one agent can efficiently track and intercept another agent, whilst 
camouflaging its own motion and minimizing its energy consumption. 

1 Introduction 
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The past decade has witnessed increasing interest in a phenomenon known as 'motion camouflage' - a 
term coined by Srinivasan and Davey (1996) to explain a class of interactive manoeuvres in hoverflies, 
observed by Collett and Land (1975). Srinivasan and Davey noted that, in certain circumstances, one 
hoverfly (A) would 'shadow' another hoverfly (B) by moving in such a way as to mimic the optic flow 
that would be produced by a stationary object in B's retina, thus 'camouflaging' its own motion and 
escaping detection. As illustrated in Fig. 1, A (the shadower) tracks B (the shadowee) by moving in 
such a way that A always remains on a straight line connecting B with a fixed, stationary point (S) in 
the environment. If A moves in this way, the trajectory of its image will be indistinguishable from the 
trajectory produced by a stationary object located at S. The lines connecting the shadowee to the fixed 
point are termed 'camouflage constraint lines' (CCLs), because successive positions of the shadower are 
constrained to lie on these lines if motion camouflage is to be maintained. This analysis assumes that 
A is small enough or sufficiently far away from B that that it appears as a point object in B's retina, 
and thus does not create any noticeable image expansion or contraction. 
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(a) Static Point Camouflage, 
static point behind 



(b) Static Point Camouflage, 
static point between 



shodowee (BJ 
(c) Infinite Point Camouflage 



Figure 1 : Depiction of the primary types of motion camouflage 



How should an animal move if it is to minimize energy consumption (for example), or interception 
time, whilst camouflaging its motion? In general, such trajectories can be calculated only if the trajec- 
tory of the shadowee is known, or can be reliably predicted [5] , [14] . Anderson and McOwen pQ tackled 
this limitation by training artificial neural networks to control simulated shadowers. Their design used 
a control system formed from three separate neural networks, and calculated a control command based 
upon the distance from the chosen static point, the direction required for camouflage maintenance, and 
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the line-of-sight error angle to the target. However the usual limitations of neural network solutions 
still applied. An exploration of the entire state-space was necessary prior to commencement of the 
learning period for the networks. The resulting controller was tailored to a specific type of shadowee 
trajectory, and lost performance quality when encountering unexpected varieties of motion. 

A more traditional approach is to create a linear quadratic feedback controller to govern the be- 
haviour of the shadowee [5]. The deviation of the shadower's trajectory from that dictated by the CCLs 
(in two dimensions) is used in conjunction with a weighted penalty on control action to create a cost 
function that results in a stable feedback controller. However, without a limit placed on the maximum 
control effort that could be requested, the required control acceleration can be very large, especially 
if a shadowee executes a trajectory that is highly non-linear. And placing a realistic cap on the con- 
trol commands compromises the accuracy of the resulting trajectory. Moreover, such a controller is 
a) not biologically realistic and b) designed to perform on a predictable shadowee path. Erratic or 
non-predictable paths are very calculation-intensive, as the solution has to be recalculated for each 
movement made by the shadowee. 

Reddy et al. explored dynamic models in both two [7] and three [13j dimensions, restricting the 
case to motion camouflage where the shadower emulates a fixed point that is infinitely far away. Their 
research was inspired by insect-capture behaviour in echo- locating bats. They modelled the shadower 
and shadowee as points subject to curvature control, rather than speed control (as is more traditional). 
They derived a feedback control law for motion camouflage using neural networks, which they called 
'motion camouflage proportional guidance' [7 , and then extended this planar guidance law to three 
dimensions, describing the particle trajectories using natural Frenet frames [7], [8]. 

Motivated by the work of [14] and [11] , Glendinning [12] developed an explicit mathematical de- 
scription of the geometry governing a motion camouflaged interaction. While an infinite number of 
potential camouflage paths exist for any given prey trajectory, his analysis allows us to prescribe further 
constraints on the behaviour of the shadower, that can give rise to a unique solution. 

Here we present a method for devising 3-dimensional trajectories which not only achieve perfect 
motion camouflage but do so in an energy-efficient manner. This approach will be useful in the design 
of energy-optimal controllers. 



2 Background 

Glendinning |12j characterized motion camouflaged paths mathematically in the following way. In 3-D 
Cartesian space, let ro(t) and rx(t) be the vectors representing the time- varying paths of the pursuer 
and the target, respectively, and let rp be the vector that represents the position of the fixed point. 
The camouflage constraint then requires that, at all times t, 

r D (t)=r D (0) + k(t)(r T (t)-r P ), 

ie 

v P -T D {t) = k{t){v P -v T {t)) V i>0 (1) 

where {fc(i)|| — oo < k(t) < 1} is a continuous time-varying scalar function. Positive k(t) corresponds 
to figure 1(a) negative k(t) corresponds to figure 1(b) (Infinite point camouflage, as in figure l(c)| 



requires a slightly different formulation, as discussed in section I3T2|) 

This equation is simply a quantitative statement of the requirement that the pursuer should always 
be located on the straight line connecting the target to the fixed point, k may be thought of as 
representing the ratio in range between the shadowee and shadower, with reference to the static point 
P, at any time t. Since k(t) can be any (bounded) continuous function, there are infinitely many 
trajectories that will satisfy the requirement of motion camouflage. The task of the analyst, then, 
is to determine constraints which lead to a unique path description. For example, we may require 
interception to occur in minimum time, or, as in this paper, with minimum expenditure of energy. 
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3 Procedure 



3.1 Camouflage with a non-infinite static point 

We assume the function k(t) to be twice differentiable, and that the initial conditions of the system 
are known. As is traditional, we ignore any gravitational effects. We then take an energy minimization 
approach to finding an optimal path, using Lagrangian mechanics. In classical and quantum mechanics, 
equations of motion are derived from what is termed the 'principle of least action'. 'Action' is a quantity 
which has dimensions of energy integrated over time. The energy in question is expressed by an equation 
known as the Lagrangian, and is generally taken as being the kinetic energy of the system minus its 
potential energy. Let the transpose of a vector x be represented as x', then assuming a particle of unit 
mass, the Lagrangian for the motion-camouflaged system can be written as 

C = KE-PE= if^fD (2) 

ie kinetic energy minus potential energy [9] (note that since we are ignoring any gravitational force, the 
potential energy term can be neglected). By solving the Euler-Lagrange equations for such a system, 
we can use Hamilton's Principle of Least Action [5] to find the path the particle will take if no other 
forces are acting on it besides those required to keep the motion camouflaged. To accomplish this, we 
set the system cost to be the integral of the Lagrangian, 



1 r* 

J=o ^Ddt (3) 



Differentiating equation ([I]) with respect to time, we obtain 

- r D (t) = fe(t)(r P - r T (i)) - k(t)r T (t) (4) 

Since k(t) parameterizes the path of the shadower, finding an appropriate stationary value for the 
energy curve 

1 

2 

is equivalent to finding the function k(t) such that the quantity 



r D (t)'r D (t)dt (5) 



to 
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k(t)r T (t)-k{t)(r P -r T (t) [k(t)r T (t) - k(t)(r P - r T (t)] dt (6) 



is minimised. We do this by solving the Euler-Lagrange (E-L) equation |5]. With C as before, the E-L 
equation for the system is 

dk dtdk ' u 
Expanding this, we obtain a second-order differential equation in k: 

k [(r P - r T )'(r P - r T )] + 2(-k [r' T (r P - r T )]) + k [-r' T (r P - r T )] = 0. (8) 

Setting (rp — r-r(i)) = a(t), we can write ([8} as: 

ka'a + 2ka'a + ka a = 0. (9) 

In other words, 



-t(ka))-a = 0, (10) 



hence we can state that k(t) will satisfy the Lagrangian conditions when the acceleration of the pursuer 
is orthogonal to the line-of-sight vector between the pursuer and target. 

This ODE does not have a general analytical solution, but it is amenable to numerical solution 
methods. Moreover, given further information about the target path, we can develop analytical so- 
lutions. In the following sections, we will primarily focus on the homogenous solution, which occurs 
when the target travels in a straight line at constant velocity, however some attention will be paid to 
other quantifiable target trajectories. 
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3.2 Camouflage at Infinity 

A similar method can be used to generate optimal paths with camouflage at infinity [14] . We again 
use a general formulation given by Glendinning [12] : 

r T -r D = k(t)e (11) 

where e is a constant vector. For the sake of simplicity we choose e to correspond with the initial 
conditions, so 

e = r T (t )-r D (t ). (12) 

Differentiating and applying the same cost function as before, we obtain the Lagrangian and develop 
a solution. Fortunately, the Euler-Lagrange equation for camouflage at infinity is solvable analytically, 
and we find the function k(t) has a general solution of the form 

ke T e = r T e + c\t + c 2 (13) 

for some constants of integration ci, c%. 



4 Examples of optimal motion camouflage equations 

With more information about the target trajectory, we can simplify the problem and discover an 
analytic solution. In this section, we present a few examples of optimal motion camouflage against 
targets with various dynamic characteristics. 



4.1 Camouflage against a target moving at constant velocity 

We begin by assuming that the target moves in a straight line at constant velocity, ie = 0. Solving 
equation (|10p is then equivalent to solving the second-order homogenous differential equation 

^(k\\a\\)=0. (14) 

Equation [14] has a solution of the form 

k\\a\\ = cxt + c 2 (15) 
for some integration constants c\ , c 2 ; therefore 

k(t) = „ Clt + C ' , . (16) 
||rp-r T (i)|| 

This represents the general solution for the optimum form of k(t) which satisfies the specified 
constraints. Clearly it encompasses the trivial solution, k(t) = V t. Explicit values for the constants 
depend on the initial conditions, so it is evident that the solution is highly sensitive to the choice of 
starting position and velocity, as is often the case with variational problems [5]. The value of c 2 is 
found as follows: assuming the initial conditions of pursuer, target and fixed point are known. We 
allow to = 0, then ko = fc(io) = fc(0) is known. By substituting into ll6[ we solve for c 2 , 

c 2 = fco(||rp-r T (0)|j). (17) 

For an engagement that ends in an interception, c\ can be similarly found from the desired terminal 
conditions, a procedure outlined later. However for an open-ended interaction, a suitable value for the 
constant c\ can still be established from the first-order initial conditions. Differentiating (|15[) gives us 

c * ■ kf {rp -Z ] (18) 



|r P -r T || ||rp-r T " 2 
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therefore by defining the initial condition fco, we obtain 



d = k (\\r P - r T (io)||) - k 



rT(to)'(rp - r T (t )) 
\\r P - r T (*o)|| 



A positive value of fco generates a trajectory with a pursuit characteristic, a negative value gives an 
escape path. Different engagements are illustrated below. Figure 2(a) shows an example involving a 
constant-velocity target. A positive value is used for k, which results in a pursuit scenario. Figure [2 (b)| 
shows a trajectory with the same starting conditions, but in which k is negative, so that the shadower 
escapes from the shadowee. 



isorj\ Shadower 



(a) Linear interaction, fc > : fco 
0.1, fc = 0.2 



(b) Linear interaction, fc < : fc > 0, fco 
0.1, fc = -0.2 



Figure 2: Optimal solutions to linear trajectories 



To demonstrate the energy efficiency of these paths, we consider two agents with the same starting 
position and velocity, camouflaging themselves against a shadowee moving along a prescribed straight 
trajectory. One agent uses the low-energy path as derived above, while the other moves in a straight 
line to the interception point (whilst still maintaining camouflage). Figure [3] shows the paths taken. 
(For simplicity, we look at only two dimensions in this instance). 
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Figure 3: Comparison between a straight interception trajectory and an energy optimal trajectory. Initial conditions, 
in cm, are as follows: X T (0) = [30, 60], V T (0) = [650, -20], X D (0) = [300, -650], Vd(0) = [9,11] 



The initial velocity vector of the straight line path was set to point in the approximate direction of 
the interception point of the optimal trajectory, as to enable a meaningful comparison. The forward 
velocity of both pursuers at the beginning of the interaction is 14 cm/s. The final speed required for 
the straight trajectory to still be camouflaged at the interception point is 1.14 x 10 3 cm/s, whereas that 
required for the optimal energy path is very much lower, 187.8 cm/s. Given these initial conditions, 
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then, to pursue a target along a straight line while remaining camouflaged, a shadower must expend 
more than 40 times the energy that would be required by a low-energy pursuit curve. 

As alluded to previously, if it is desired to intercept the shadower within a predetermined time (a 
so-called 'finite horizon' capture), or indeed if another geometric condition is to be met, then we can 
determine the necessary constant C\ from the end boundary conditions. In the case of capture in some 
finite time tf, we say k(tf) = 1, so substituting this condition and equation 1171 into equation 1161 we 
can write: 

ci = r [llrp - r T (*/)ll ~ Mrp - M*o)||] • ( 19 ) 
Figure S] demonstrates a finite horizon capture path. 



ST 



Figure 4: Capture trajectory for a shadower with initial conditions fco = 0.1, fco — 0.2, static point at rp = 
[200,-650,500]. Shadowee initial conditions r T (0) = [30,60,150], constant velocity r T = [200,-20,60]. Capture 
horizon is 12s (CCLs shown here every 0.4s). 



A similar method can be used to find a trajectory for camouflage at infinity against a target moving 
at constant velocity. Recall equation 1131 We can again find the constants c\ and c-i from the initial 
position and velocity conditions. Alternatively, for a capture trajectory, we can determine the optimal 
Ci for a given final time from the terminal boundary conditions, so 

ci = T {Ht f )e T e - $(t f )e - c 2 ) . (20) 

For a capture trajectory, k(tf) — 0, and this is demonstrated in Figure[5] For tracking at a constant 
distance d = \\rx — fd II = ll e ll' we se * ^(*o) = k(t) = k(tf) = 1. An example of the latter tracking 
result can be seen in Figure [SI 



Figure 5: Capture using camouflage at infinity. Initial conditions and velocity of the shadowee are rxito) 
[-30150150], v T (*o) = [200, -20, 60]. Starting point of shadower is at [0, 0, 0]. 
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Figure 6: Tracking at a constant distance. Initial conditions and velocity of the shadowee are identical to Figured 
Starting point of shadower is at [0, 0, 0]. 



4.2 Camouflaging against a target following a quasi-3D guidance law 

In conspecific two-body guidance, it may be surmised that both participants have similar capabilities 
and are using similar guidance techniques. In such a scenario, it is unlikely that one participant will 
maintain a constant linear (or angular) velocity while the other does not, so the conditions on target 
dynamics assumed in the previous examples may not be useful. However in a conspecific situation, we 
may be able to use the similarity of the participants to come up with some other constraint on the 
target motion. One such example follows here. 

Recall that our main result in equation \W\ shows that the pursuer's acceleration under Lagrangian 
camouflage conditions will be orthogonal to the line-of-sight vector between pursuer and target. An 
interaction of this nature will produce a relative motion vector that lies in a plane. Bodies with tra- 
jectories satisfying this condition are said to be moving in a 'quasi-3dimensional' manner |15j . Certain 
commonly used guidance laws which result in quasi-3D motion have as their defining characteristic 
a relative acceleration of zero along the radial axis between shadowee and shadower, as does motion 
camouflage. We examine a situation where the shadower, D, is following a motion camouflaged path, 
and the shadowee, T, is following a quasi-3D guidance law - in other words, the shadower and shad- 
owee are using similar (but possibly non-identical) guidance strategies, both accelerating in a direction 
orthogonal to the line-of-sight. Using this information, we can make a general statement about the 
function k{t) under these conditions. 
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Figure 7: Geometry of motion camouflage in three dimensions in a spherical reference frame, with origin at the 
static point. 



We examine a camouflaged engagement in a spherical reference frame (Figure [7]). The origin is 
fixed at the static point of the engagement, and the co-ordinates of the pursuer and target can be 
written (rjj, 8jj, <Pd) and {tt,9t,4>t) respectively. It is straightforward to show that under motion 
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camouflage conditions, 6d = &t = 0, 4>d = <t>T — <fi- Hence we can define an orthogonal set of unit 
vectors (e r ,ee,e0) along the co-ordinate axes, as in Figure [71 

The displacement and velocity of the pursuer in this reference frame can be written 

r D = r D e r 

r D = f D e r +r D dcos(j>e e +r D 4>e (j> (21) 

hence the acceleration components along the co-ordinate axes are 

a-D r — — Td4> — r^O 2 cos 2 4> (22a) 

CLD g — rodeos 4> + 2rjj6 cos (f> — 2rr>4>6 sin <j> (22b) 

ajj^ = ro<j> + 2r d4> + rod 2 cos <j> sine/). (22c) 

According to equation ([TO)) , the acceleration of the pursuer D is entirely orthogonal to the line-of-sight 
vector e r , or in other words, the radial component of the acceleration is zero. Hence we can write 

r D - r D (0 2 cos 2 + 4> 2 ) = 0. 

Applying equation (TT]), we establish the following condition for k(t): 

^-(kr T ) = r T (e 2 cos 2 ^ + <i) 2 ). (23) 
The angular velocity of D is found in the following way: 

f2 = rD 2 V ° = —(f>eg + #cos0e^, (24) 

r D 

hence we can replace (f2"3")l with the following general equation: 

d 2 

— (kr T ) = r T n 2 . (25) 

Let r = yt — yd- If the shadowee is using a quasi-3D guidance law to react to the shadower, then 
the following holds true [T5] : 

P = m 2 (26) 

We can write the left-hand side of this equation in terms of k: 

— r'x — (krx + 2kfx + kfx) 

Doing the same to the right-hand side, we find 

^(rr) ~ ^(kr T ) = r T n 2 - (fcr T )fi 2 . (28) 

Recall the general equation described in (|2l)]) . We can rearrange this to find a differential equation in 
k, rr to describe fl: 

n2 = itS {krT) (29) 

Substituting into (|28p and expanding, we find that k satisfies the following differential equation: 

kr T + 2kr T = 

which has a solution 

r* i 

k{t) = Cl +c 2 —r^dt. (30) 

J to r T(t) 2 

The exact form of k(t) and the constants Ci, ci depend on the precise formulation of the guidance 
law being used by the shadowee and the initial conditions of the engagement. An example of two-body 
guidance using a specific quasi-3D guidance law is given in section 15.21 
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5 Possible Applications 



In order to demonstrate a possible use for these optimal path equations, we consider what kind of 
guidance system may be easily modified to produce trajectories which satisfy the above conditions. A 
common tracking strategy which has been observed in insects such as blowflies and houseflies [4] , [10] , 
is for the pursuer to apply an acceleration proportional to the relative angular velocity between the 
pursuer and target. This type of guidance is conventionally known as proportional navigation (PN), 
and as well as being observed in nature, has been used for many years in autonomous and independently 
guided tracking systems. Motion camouflage is a natural candidate for a modified PN guidance law 
[7], since the relative angular position is a key variable when it comes to defining and controlling for 
motion camouflage. 

Consider the problem of static point motion camouflage in two dimensions. We use a polar reference 
frame, as in section 14.21 however for simplicity we will restrict ourselves to two dimensions. This co- 
ordinate system can be described by (e r ,e#), where e r is a unit vector along the line-of-sight (LOS) 
between agents D and T, and eg is a unit vector orthogonal to e r . 

Given equations ([26]) and (fT0|) we can write the acceleration of the shadower D as 

a D = a De = (kr T )0 + 2^(kr T )e. 

It is straightforward to see that the above is equivalent to 

a D = ka T + 2kr T 0. (31) 

This is a form of augmented PN guidance, which we term MCPN, and depending on the value of <zt, 
this augmented equation may reduce further to a simple PN function. For example, if the shadowee is 
moving with a constant velocity, we can write equation (|3 1 1) in terms of the relative distance and line 
of sight velocity, 

kr ■ 

fto = 2 



1-k 

= TrO (32) 

where T is a variable gain function depending on the ratio (k) of the distances of the shadower and 
shadowee from the static point. 



5.1 MCPN Guidance against a target moving at constant velocity 

We assume an oblivious target moving at constant velocity. The relevant guidance law is described in 
equation [32] Figure [8(b)] shows the camouflage path achieved using a PN-guided pursuer following a 
motion camouflaged path under locally energy- minimal conditions (MCPN) , and the radial and angular 
acceleration required. 



5.2 MCPN guidance against a reactive guided target 

Previously ( £|4.2p we discussed a conspecific pursuit scenario. There, we derived a motion camouflage 
equation that may be used against a reactive shadowee which is simultaneously tracking the shadower. 
This may possibly be of use to robotics researchers or autonomous agents - it will almost certainly 
give us some insight into the flight patterns seen in insects known to use motion camouflage. Hence 
we now develop a guidance law that may be used to maintain camouflage against a similarly capable 
manoeuvering target, provided that some rudimentary facts about the target dynamics are known. 

Suppose the shadowee T has a non-constant velocity, and in fact is also tracking the pursuer, using 
a known guidance law. Let the shadowee be using a quasi-3D proportional guidance law known as True 
Proportional Navigation (TPN), which takes the form 

a T = Xr 9. (33) 
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(a) Shadower and shadowee trajectory, 
given initial conditions: ko — 0.2, ko = 0.1, 
r T (0) = (30,60), r T (0) = (200,-20), r P = 
(200, -650) 



(b) Acceleration commands in the spherical 
frame (A r e r ,A t ee) 



Figure 8: Motion Camouflage path for a non-maneuvering target using PN-derived acceleration commands. Note 
that in figure (b), A r = 0, as expected. 



In this formulation, A is a gain constant, tq is the initial range between the shadower and shadowee, 
and 9 is the line-of-sight vector between shadower and shadowee (from the shadowee's perspective). 

Recall that according to Yang and Yang [15], with a quasi-3D guidance law, the relative radial 
acceleration an r — 0. From (|10| we know an r = 0, hence we can state that ax r = 0. Therefore for a 
shadowee using TPN against a motion camouflaged shadower, the following two equations hold: 

r T = r T 9 2 (34a) 
a T = r T 9 + 2f T 9 = \r 6 (34b) 

We make the following variable substitution [TS]: let V r = ry, Vg = ttO, then the equations in |3"4"1) 
become 

V T - Vg9 = (35a) 
V e + V r 6 = Xr Q 9 (35b) 

We can change the independent variable t to 9 and solve the above to find 

V r = Acos(9 + B) + Xr a (36) 
Vg = -Asm(9 + B) (37) 

where A, B are constants (which can be found from initial conditions). From (|3~T|) and (|3T)|) . k(t) takes 
the following form: 

(38) 




cos(6» n + B) 

which can be solved numerically. The shadower acceleration can be written 



(39) 



a D = (Xkr + 2rj^)9 = T*9, (40) 
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which again gives us a PN guidance law, this time with variable gain T*. A guided shadower using a 
low-energy path derived with the above method against a guided target can be seen in Figure 9(a)[ 
and the desired accelerations of both target and pursuer are shown in Figure |9(b)| 



200 400 



1000 1200 1400 1600 1600 2000 




(a) The shadowee path is labelled 'TPN 
law', the shadower's path is labelled 
'MCPN'. Initial conditions are: rp — 
(2000 - 1650), r T (0) = (30,60), r T (0) = 
(200, -20) . 



(b) Acceleration commands in the respec- 
tive body-centred spherical frame for both 
shadowee (TPN) and shadower (MCPN) 



Figure 9: MCPN path against an active target using TPN guidance, (a) Pursuer and target trajectory, given 
ka = 0.08, fco = 0.1 (b) Acceleration commands in the spherical frame for both pursuer and target (A r e r , A t eg). 
For both bodies, the radial acceleration is zero. As the endgame approaches, the acceleration of both shadowee and 
shadower increases, but over the course of the interaction the acceleration of the MCPN body is notably less than 
that of the TPN body. 



It is hoped that by finding elegant and efficient means of mimicking pursuits and interactions 
observed in nature, we can gain insight into the guidance mechanisms used by insects and other 
animals. For example, research by Boeddeker, et al [I], and Land, et al [TO], has found that insect 
pursuit systems can be modelled quite simply using a static first-order gain controller with only a few 
input variables readily available via the visual system. Similarly, Ghose, et al [B], show evidence that 
bats use an interception strategy commonly found in guided missiles to intercept prey. By combining 
these simple models with a guidance system that enables us to replicate gross characteristics of known 
flights, we can obtain a streamlined but fully controllable navigational system. 



6 Conclusion 

This study describes how energy-efficient strategies can be derived for tracking or shadowing moving 
objects whilst maintaining motion camouflage. A series of control strategies is presented for implement- 
ing energy-efficient motion camouflage under various conditions. The findings can be usefully applied 
to a variety of tracking situations, as well as provide insights into some of the behaviours exhibited by 
insects and animals. Our most important result was the discovery of a necessary orthogonal condition 
to the pursuer acceleration under an energy-optimal motion camouflage path. 

This work was funded by an ANU Postgraduate Scholarship (to N.C.) and had partial support from 
the ARC Centre of Excellence in Vision Science (to M.S.). The authors would like to thank Jochen 
Zeil for discussion and advice and Alan Carey for an illuminating conversation. 
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